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1.   Summary :   It  sometimes  happens  that  a  physical  phenomenon  (target) 
leaves  occasional  traces  of  its  existence  as  patches  of  emitted  material 
(blobs)  that  may  temporarily  grow  in  size,  but  which  ultimately  dissipate 
and  become  undetectable.   In  this  report  we  investigate  the  chances  of 
contacting  one  of  these  blobs,  assuming  that  the  target  moves  in  its  medium 
enough  to  prevent  effective  overlap  of  one  blob  on  another.   Search  is  assumed 
to  be  random  within  some  confined  region  by  a  sensor  capable  of  detecting  a 
given  concentration  of  emitted  material;  the  effect  of  this  assumption  is  to 
make  the  contact  rate  at  any  time  proportional  to  the  total  projected  length 
(in  two  dimensions)  or  area  (in  three  dimensions)  of  all  extant  blobs.   The 
random  search  assumption  is  made  for  analytical  convenience;  nonetheless, 
the  assumption  is  often  a  surprisingly  accurate  representation  of  real 
searches,  particularly  when  the  target  is  moving. 

Let  C  be  the  mean  number  of  times  that  a  blob  is  contacted.   Formulas 
for  C  are  derived  in  sections  3  and  4  for  the  cases  of  diffusion  in  two  and 
three  dimensions.   As  long  as   C  is  much  smaller  than  1,  the  mean  time  E(T) 
until  the  first  blob  is  contacted  is  approximately  an  exponential  random 
variable  with  mean  1/XC,  where  X     is  the  Poisson  rate  at  which  blobs  are 
emitted.   Actually,  1/XC  is  a  slight  underestimate  of   E(T),  as  is  shown 
in  sec.  2.   Exact  formulas  for  the  probability  distribution  of  T   (formula 
A4)  and  E(T)   (formula  A5)  valid  for  all   C  are  derived  in  the  appendix. 

The  blob  that  is  detected  may  be  an  old  one  if  blobs  have  a  long 
lifetime.   The  time  required  to  detect  a  blob  whose  age  does  not  exceed   t 
could  be  analyzed  by  truncating  the  rate  of  detection  function  at 

t  •   i  e   r(x)  =  0   for  x  >  t  .   The  effect  would  be  to  reduce  the 
o'    "  * '  o 

dimensionless  constant  in  formulas  (9)  and  (16) . 


2 .   The  mean  time  to  detection  exceeds  1/XC 

On  account  of  the  random  search  assumption,  the  total  number  of  contacts 
when  the  blob  emission  times  are  given  is  a  non-homogeneous  Poisson  process. 
Let  Y(t)   be  the  mean  value  function  of  this  process,  and  let  X(t)   be 
the  actual  number  of  contacts  up  to  time  t.   Since  Y(t)   is  itself  a 
stochastic  process,   X(t)   is  doubly  stochastic  [1].   Since  X(t)   is  a 
Poisson  random  variable  with  mean  Y(t),   E(x(t)|Y(t))  =  Y(t),  and 
E(x(t)J  =  E(Y(t)J   by  conditional  expectations.   This  is  handy  because 
E(Y(t)J   is  easy  to  predict  in  equilibrium.   In  fact,  let  C  be  the  average 
number  of  contacts  per  blob  (see  sections  3  and  4  for  examples  of  how  C 
can  be  computed).   Assuming  that  blobs  are  emitted  at  the  Poisson  rate  X, 
we  must  have 

(1)  E(x(t))  =  E(Y(t))  =  XCt 

It  follows  that  the  average  time  between  contacts  is   1/AC.   However,  this 
is  not  the  same  thing  as  the  mean  time  E(T)   until  the  first  contact, 
since  the  origin  of  T  is  meant  to  be  an  arbitrary  point  in  the  process, 
rather  than  a  point  at  which  a  detection  has  just  occurred  and  for  which 
Y(t)   is  therefore  likely  to  be  large. 

00 

Since  T  is  a  positive  random  variable,   E(T)  =  /  P(T  >  t)dt,   and 

oo  0 

therefore  E(T)  =  E(/  P(T  >  t|Y(t))dt)   by  conditional  expectations.   But 

0 
the  event  (T  >  t)   is  the  same  as  the  event   (X(t)  =  0) ,   and 

P(X(t)  =  0|Y(t))  =  exp(-Y(t)),  so 

f  °°  \    °° 

(2)  E(T)  =  E  /  exp(-Y(t))dt   >  /  exp(-E(Y(t) )  )dt 

^0  >         0 

00 

=  /  exp(-ACt)dt  =  1/AC  , 
0 


where  the  inequality  is  Jensen's  inequality  applied  to  the  convex  function 
exp(-x).   So  the  mean  time  to  the  first  contact  exceeds   1/XC,  in  general. 
The  condition  for  E(T)   to  be  approximately  1/XC  is  that  Y(t)   should 
not  be  "very  random".   Intuitively,  this  will  be  the  case  if  the  contribution 
due  to  a  single  blob  is  small;  i.e.,  if  C  <<  1.   This  is  correct,  but  the 
proof  is  not  trivial,  so  we  defer  it  to  the  appendix. 

3.   Evaluation  of  C   for  turbulent  diffusion  in  the  ocean 

In  [A],  it  is  stated  that  the  diffusion  of  substance  from  a  point  in  the 

ocean  is  essentially  a  two  dimensional  problem,  since  there  is  very  little 

vertical  mixing  compared  to  turbulent  horizontal  processes.  Several  radially 

symmetric  formulas  are  offered  to  describe  the  way  in  which  the  concentration 

of  a  substance  s(t,r)  changes  with  distance   (r)   and  time  (t)   from  the 
point  of  injection,  most  of  which  are  special  cases  of 

(3)  s(t,r)  =  ^jr exp(-rb/(St)a), 

TT(3t)ZCt/br(l+2/b) 

where  M  is  the  amount  of  material  released,   a  and  b  are  dimensionless, 
and   $  has  whatever  units  are  required  to  make   s(t,r)   have  units  of 
material  per  unit  area.     T  is  the  Gamma  function,  and  is  required  in  (3) 
in  order  to  make  the  volume  underneath  s(t,r)   equal  to  M  for  all   t. 
The  case  of  classical  diffusion  is  a  =  1,  b  =  2. 

At  time   t  after  injection,  a  concentration  sensor  will  be  able  to 
detect  the  substance  within  some  distance  p(t)   of  the  point  at  which  it  was 
released.   Suppose  that  M  units  of  material  could  just  barely  be  detected 
inside  a  circle  of  radius  r_   if  the  material  were  spread  uniformly;   r 

is  a  measure  of  the  sensor's  sensitivity.   Then  p (t)   must  satisfy 

2 
s(t,p(t))  =  M/(7Tr  )   at  all  times   t.   Using  (3)  to  solve  for   p(t),  we 

obtain 


(4) 


P(t)  = 


(3t)a/b  £n 


(et)2a/br(i+2/b) 


1/b 


except  that  p(t)   is  taken  to  be  0  rather  than  negative.   If  the  sensor 
is  conducting  a  random  search  at  speed  v  in  an  area  of  size  A,  then  the 
rate  of  contacting  the  blob  at  time  t  is  2vp(t)/A,  since  2p(t)   is 
the  projection  of  the  circular  blob  normal  to  the  searcher's  track.   The 
average  number  of  detections  per  blob  is  therefore  ♦ 


(5) 


C  =  /  (2vp(t)/A)dt  , 
0 


Let 


K  =  r^/rd+2/b), 


2a/b 


substitute  p(t)   from  (4)  in  (5),  and  then  substitute  u  =  (St)    /K. 
The  result  is 


(6)     C  =  (2v/A)  /  [Ku]1/2[£n  ^]1/b 

0 


Kb  du 


2a3(Ku) 


l-b/2a 


where  the  upper  limit  of  u  is  the  largest  number  for  which  the  integrand  is 
positive.   Factoring  out  constants  and  collecting  powers  of  Ku ,  (6)  is 


(7)     C  -  2&  /  [KU,1^-1/2  [to  1,1/*  du 


Aa8 


0 


The  definite  integral  has  a  closed  form  expression  [  3]  .   Substituting  this  in 
(7) ,  we  obtain 

(8)     C  =  vKb/2a+1/2  b   ra+l/b)/[AaB  (b/2a+l/2)1+1/b]  . 


Recalling  the  definition  of  K,  we  finally  have 


v  rb/a+1  bru+l/b) 
(9)  C  =  


AaB  r(l+2/b)b/2a+1/2  (b/2a+l/2)1+1/b  ' 

In  classical  diffusion,  since  T(3/2)  =  /F/2  and  T(2)  =  1,  (9)  re- 
duces to 

v  r  fn 

(10)  c  = ^—TP)  - 

AB  (3/2)J/Z 

However,  classical  diffusion  is  not  a  good  description  of  turbulent 
spreading  in  the  ocean.  According  to  reference  [4],  better  descriptions 

would  be  a=2,  b=2,  8=2  cm. /sec.  (Okubo  and  Pritchard) ,  or  a=3,  b=2, 

2/3 
B=.05  cm   /sec,  (Okubo).   Both  of  these  models  are  from  table  5  of 

reference  [4].   If  v=10  cm/sec,  r  =10  cm.,  and  A=10  cm  ,  (9)  gives 
C  =  .44  for  the  first  of  these  models,  and  C  =  1.55  for  the  second.   These 
calculations  may  give  the  reader  some  idea  of  circumstances  where  C  «  1. 
The  fact  that  different  investigators  have  arrived  at  models  that  provide 
answers  differing  by  a  factor  of  3.5  should  not  be  disturbing  when  the  sub- 
ject is  turbulence  in  the  ocean. 

Figure  1  shows  (4)  for  the  same  two  models.   The  maximum  blob  size 
depends  only  on  r   and  b,  and  is  consequently  the  same  in  both  cases.   The 
fact  that  the  area  under  one  curve  is  3.5  times  as  large  as  the  area  under 
the  other  is  essentially  a  question  of  blob  lifetime. 

4.   Classical  diffusion  in  three  dimensions. 

This  derivation  will  parallel  that  of  sec.  3.  If  S(t,r)  is  now  the 
concentration  of  material  per  unit  volume  and  D  the  diffusion  constant, 
then  (3)  is  replaced  by 

(11)  S(t,r)  = M  -/9   exp(-r2/2Dt). 

(2TTDt)J/ 

We  also  have 

5 


(12) 


(4/3)7rr 

-£p>   exp(-p2(t)/2Dt), 


(2-rrDt) 


3/2 


where  rn   is  the  radius  of  the  largest  sphere  within  which  M  units  of 

material  would  be  detectable  if  spread  evenly. 

2  1/3 

Solving  (12)  for   p  (t),  and  defining  K  =  (9tt/2)    ,   we  have 


(13) 


p  (t)  =  2Dt  £n 


2  3/2 


KDt 


=  3Dt  in 


0 


KDt  ' 


except  that   p(t)  =  0   for   t  >  r n/KD.   The  rate  of  detection  at  time   t 

2        "\ 

is   now     r(t)   =  v[Trp    (t)J/V,      where     V     is    the  volume  of   the  confining  region, 

so 


(14) 


C   =   ttv/V   /      p    (t)    dt. 
0 


Substituting  u  =  KDt/r  , 


(15) 


C  =  (Tiv/V)(r^/KD)(3r^/K)  /  u  iin(l/u)  du 
u       u    0 


The  definite  integral  [3]  is  1/4,  so 


(16) 


C  =  .1343  (vr^/VD) 


Appendix: 

Let  r(x)   be  the  detection  rate  on  a  blob  at  a  time  x  after  it  has 

been  emitted,  and  suppose  blobs  are  emitted  in  a  Poisson  process  with  rate  A. 

Let  N(t)   be  the  number  of  blobs  emitted  in  [0,t],  and  let   S.   be  the  times 

at  which  they  are  emitted,  i  =  1,  ...,  N(t).   For  0  _<  u  _<  t,  the  total 

N(t) 
detection  rate  at  time  u  is   £   r(u-S.),  where  r(x)   is  taken  to  be  0 

i=l       X 

for  x  <  0.   Therefore  the  average  number  of  detections  in  the  interval 

[t-A,t],  given  the  emission  times,  is 

N(t)   t 
(Al)  Y(t)  =  I        /  r(u-S.)  du,   A  _>  0. 

i=l  t-A      x 

The  stochastic  process  Y(t)   is  a  filtered  Poisson  process  with  impulse 

t  t-T 

response  h(t,x)  =  /  r(u-T)  du  =  /    r(v)  dv.   The  generating  function 

t-A  t-T-A 

of  Y(t)   is  therefore  known  [2]: 


(A2)  ^(S)  =   E(exp(SY(t))  =  exp 


t 
X  /   [exp(Sh(t,T))  -  1]  dT 
0 


By  substituting  x  =  t  -  T,  this  can  be  written 


(A3)  <J/t(S)  =  exp 


t         x 
A  /   [exp(S  /  r(v)  dv)  -1]  dx 
0        x-A 


The  number  of  detections  in  the  interval   [t-A,t]   is  a  Poisson 
random  variable  with  mean  Y(t).   The  probability  of  no  detections  in  this 
interval  is  therefore  E(exp(-Y(t) ) )  =  i(>  (-1) .   If  T  is  the  time  from 
t-A  to  the  next  detection,  we  therefore  have  P(T  >  A)  =  \p    (-1).   Letting 


C  •*■  °°  to  obtain  equillibrium, 


(A4) 


P(T  >  A)  =  exp 


X  /   [exp(-  /  r(v)  dv)  -1]  dx 

0        x-A 


Since  T  is  a  positive  random  variable,  we  also  have 

00 

(A5)  E(T)  -  /  P(T  >  A)  dA, 

0 

and  this  provides  an  explicit  formula  for  E(T)   that  could  be  evaluated  by 

numerical  integration  in  any  specific  case.   However,  suppose  the  mean 

00 

number  of  contacts  per  blob  (C)   is  «1.   Since  C  =  /  r(v)  dv,  it  follows 

0 
x 

that   /  r(v)dv  «1,  and  therefore  that 

x-A 


1  -  exp 


x        1 

-  /  r(v)  dv  %  /  r(v)  dv,   and  therefore  that 

■  x-A       J   x-A 


(A6) 


P(T  >  A)  %  exp(-X  f(A)),    where 


(A7) 


f(A)  =  /   [  /  r(v)  dv]  dx 

0   x-A 


Since   f (0)  =  0,  and  since   (d/dA)  f(A)  =  /  r(x-A)  dx  =  C,   it  must  be  true 

0 

that   f(A)  =  CA,  and  therefore  that   T  is  approximately  exponential  with 
mean  1/XC.   Furthermore  since  all  of  the  above  approximations  can  be 
written  as  inequalities, 


(A8) 


E(T)  >   1/XC, 


with  equality  holding  when  C  <<  1, 
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